Andreev Bound States in the Kondo Quantum Dots Coupled to Superconducting 

Leads 
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We have studied the Kondo quantum dot coupled to two superconducting leads and investigated 
the subgap Andreev states using the NRG method. Contrary to the recent NCA results [Clerk and 
Ambegaokar, Phys. Rev. B 61, 9109 (2000); Sellier et al., Phys. Rev. B 72, 174502 (2005)], we 
observe Andreev states both below and above the Fermi level. 
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When a localized spin (in an impurity or a quantum 
dot) coupled to BCS-type s-wave superconductors(S), 
two strong correlation effects compete. The supercon- 
ductivity tends to keep the conduction electrons in singlet 
pairs0|, leaving the local spin unscreened. The spin state 
of the total wave function is thus a doublet. The Kondo 
effect tends to screen the local spin with the spins of the 
quasiparticle excitations in the superconductor. The to- 
tal spin state is thus a singlet 0|- This competition gives 
rise to a quantum phase transition from the doublet to 
singlet state, and the transport properties are dramati- 
cally changed across this transition^ 0. For example, 
the Josephson current through the quantum dot (QD) 
coupled to two superconducting leads has a 7r-shift in its 
current-phase relation (a so-called 7r-junction behavior) 
for the doublet state, while it has a 0-iunction behavior 
for the singlet state |1 II II 1111 EJuUliJ^^ fact, 
all the previous studiesllllBB HMinnHlli of the 
transition between the doublet and singlet state focused 
on the current-phase relation Is {(p) of the Josephson cur- 
rent. 

However, there are another non-trivial issues, about 
the Andreev bound states in such a system. The issues 
are: (1) How many subgap Andreev states are there? 
(2) Are the subgap Andreev states true bound states or 
quasi-bound states with finite level broadening? 

Using the non-crossing approximation (NCA), Clerk 
and Ambegaokar 7] investigated the close relation be- 
tween the O-TT transition in Isifj}) and the Andreev states. 
They found that there is only one subgap Andreev state 
and that the Andreev state is located below (above) the 
Fermi energy Ep in the doublet (singlet) state. They 
provided an intuitively appealing interpretation that in 
the doublet state, the impurity level well below Ep is 
singly occupied and due to the strong on-site interac- 
tion energy U only hole-like excitations are allowed. In 
the singlet state, a small probability to find the impu- 
rity empty, electron-like excitations are allowed. This 
result was supported further by a more elaborated NCA 
method by Sellier et al. As stressed by Clerk and 

Ambegaokar 0, this observation casts a strong contrast 
with the non-interacting case, where bound states always 
occur in pairs (below and above -Ef)E3,0|- Moreover, 



the Andreev state has a finite broadening. 

To the contrary, the Hartree-Fock theory (HFT) [3| 
predicts two Andreev states both below and above Ep 
in the Kondo regime. This was in agreement with the 
numerical renormalization group (NRG) calculations by 
[l7j , who studied the Andreev states as a function of the 
impurity level position. Slave-boson mean-field theory 
(SBMFT) [i3 also predicts both electron-like (above Ep) 
and hole-like (below Ep) Andreev states. Further, they 
both predict infinitely sharp Andreev states. However, 
HFT and SBMFT are effectively non-interacting theory 
and may not provide a strong argument against the NCA 
result. In a previous workjgj we observed both Andreev 
states in the NRG result, but the model had the particle- 
hole symmetry. 

In this work, we report a systematic study of the is- 
sues using the NRG method. We find both electron-like 
and hole- like Andreev states except in the deep inside the 
doublet phase (superconducting gap even bigger than the 
hybridization) . We also provide a supporting interpreta- 
tion based on the variational wave functions suggested 
by Rozhkov and Arovas |^. 

The subgap Andreev bound states are important be- 
cause they are directly related to the transport properties 
of the S/QD/S systems as in the recent experiment 0, 
mi^lill' Recent developments in mesoscopic transport 
experiments may allow direct measurement of these An- 
dreev states through tunneling spectroscopy. So far, the 
NCA/SNCA and the NRG are the only methods which 
can treat rather systematically the many-body correla- 
tions of the Anderson-type impurity coupled to super- 
conductors. The discrepancy between the two methods 
may motivate further theoretical efforts for better under- 
standing of the many-body states of the system. 

Summary of the Results We consider a QD (or mag- 
netic impurity) coupled to two superconducting leads. 
The Hamiltonian H — He + Hp, -I- Hp consists of three 
parts: Hq describes two, left (L) and right (i?), BCS- 
like s-wave superconductors with the superconducting 
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gap ^L{R) and band width D 



Hc = 



E 

l=L,R L ka 



(1) 

We assume identical superconductors, £,Lk — C-Rfe — & 
and = = Ae+"'^/^, where is the phase differ- 
ence between the two superconductors. iJ/j describes an 
Anderson-type single level in the QD 



(2) 



with Tier = dl-da- f-d is the single-particle energy of the 
level and U is the on-site interaction. Here we consider 
particle- hole asymmetric case (e^ ^ —U /2). To compare 
our results directly with the NCA results, we will mostly 
take L/ = oo by preventing double occupancy. Finally Ht 
is responsible for the the tunneling of electrons between 
the superconductors and the QD 
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(3) 



For simplicity we will assume the symmetric junctions 
with tunneling elements insensitive to the energy, V^fc = 
Viik ~ V . The broadening of the level is given by F = 
ttPl{Ef)\Vl? +t^Pr{Ef)\Vr\^ = 2np{EF)\V\^ . Yov 
calculation, we followed the standard NRG method 23, 
MUEHEIil extended to superconducting leads [l7 ll2a 

ig. 

There are two competing energy scales in the system. 
The superconductivity is naturally governed by the gap 
A. The Kondo effect is characterized by the Kondo tem- 
perature Tk, given by ^8, 11, 30, 31] 
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where W^o = min{Z), [/}. For Tk ^ A, the ground state 
is expected to be a singlet and the Joscphson current is 
governed by the Kondo physics. In the opposite limit 
Tk <C a, the ground state is a doublet and the trans- 
port can be understood perturbatively in the spirit of 
the Coulomb blockade (CB) effec 
happens at Tk ~ A. See Fig. |2| 

Figure ^ summarizes the results. Figure ^ (a) shows 
the positions, E^, and E^, of the subgap Andreev states 
ioi U — oo and = 0. We observe two Andreev states, 
below and above Ep are observed in a wide range of 
A/Tk (in particular on both sides of the transition point 
Ac/Tr- ~ 1). More important are Fig.n](b) and (c), the 
spectral weights Ae [At) of the electron- like (hole-like) 
Andreev states, defined by 
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FIG. 1: (color on-line) The energies (a) and the corresponding 
spectral weights (b) of the subgap Andreev states for U — oo 
Anderson model, (c) The ratio Ae/Ah of the spectral weights 
of the Andreev states. Inset: Ae/Ah as a function of —U /2ed- 
(d) Average occupation of the QD level, ed — —O.ID, U — oo, 
r = 0.02D {Tk ^ 3.88 x 10"^). 



near E ~ Ep {p — e,h). Except for very large A (A ^ 
F), both Ee and Eh have order of magnitude the same 
spectral weights. These observations are consistent with 
the occupation of the dot level shown in Fig.n](d). Unlike 
the intuitive interpretation by the occupation does 
not change much across the transition point, although 
there is a small jump [emphasized in the blue circle in 
Fig. n (d)] . The results in Fig. ^ remain qualitatively 
the same for finite U; an example is show in the inset 
of Fig. ^(c). Finite phase difference (not shown in the 
figure) does not make any qualitative change (about the 
existance of the Andereev states both above and below 
Ep), either. 

Variational Theory The main features of the results 
summarized above can be understood qualitatively in 
terms of the variational wave functions [H . For the singlet 
state in the U = oo limit we take the trial function of the 
form 

I ^"^qeL^R 

+ j:C..'i,ll-]\0) (6) 



with Cqq' = Cq'q. FoT the doublct state we take 

\D,) = \Ad\ + Y.B,il, + T.^'i'i'iii'A 

I q qq' 

)^T.^'^'''ii{i'^d,-jl,,dM\o) (7) 
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FIG. 2: (color on-line) Results from the variational calcula- 
tions. Plotted is the Ratio Ae/Ah of the spectral weights as 
a function of — ed/F for various values of F. A/D = 0.1. 

with Cqqi = Cqiq and Dqqi = —Dqiq, and analogously 
\Di). The coefficients in the trial wave functions are 
determined by minimizing E = (*|7Jvl/) /(^'|vl') for \^) = 
\S) and I*) = \D,). 

From the form of the trial wave functions, the spec- 
tral strength of the hole-like Andreev state in the singlet 
phase depends on how much the QD is occupied (Bq) in 
\S) and how much the QD is empty (Bq) in \Da-), namely, 
on the matrix element 



(up to normalization constant y/ {S\S) (Z)| \D^)). The 
strength of the electron-like Andreev state in the doublet 
phase also depends on {Di\di\S). Likewise, the strength 
of the electron-like (hole-like) Andreev state in the sin- 
glet (doublet) phase is determined by the matrix element 

{D^\d]\S)=A*A + Y,C*qq,Cqq,. (9) 
qq' 

According to the NCA resultsjl E^, {D^\di\S) 
((D| I (ills')) should vanish in the singlet (doublet) phase. 
However, as shown in Fig. |21 neither of them vanishes, 
and the spectral weights Ag and Ah are similar in or- 
der of magnitude on both sides of the transition point, 
in agreement with the NRG results. We must point out 
that the agreement between the variational and NRG re- 
sults is only at a qualitative level. The ratio Ae/Ah from 
the variational method is about 5 times bigger than the 
NRG result. However, this is not surprising because the 
variational method is limited in the regime Tr- not too 
large compared with A; see below. 

There is another interesting point to be noticed in 
the variational wave functions in Eqs. © and 0. The 
lowest-energy solution to the variational equation for \S) 
is well separated from the continuum. This is also true 
for \Da-)- This suggests that the subgap Andreev state 



FIG. 3: (color on-line) Phase diagram of the Anderson impu- 
rity model with superconducting leads (for (p = 0). The phase 
boundaries for the infinite-?7 model (red solid line with cir- 
cles) and for the particle-hole symmetric model (black dashed 
line with squares), respectively, have been calculated by the 
NRG method. 
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FIG. 4: (color on-line) Spectral weights A^ (^h) of the 
electron- like (hole-like) Andreev state as a function of A/Tk 
at = = 0. Inset: the same in a wider range of A/Tk 
including the doublet phase. 

is a true bound state without broadening. Wc will come 
back to this discussion below. 

Universality Since the singlet-doublet transition in 
the system is a true quantum phase transition, the uni- 
versality is also an important issue. With A and 
being the only two low energy scales in the system, phys- 
ical quantities should depend only on the ratio of A/T^ 
but not on the details of the system. 

In Fig. 0] we plotted the normalized spectral weights 
Ae{h) /A as a function oi A/Tk for various values of e^/r. 
We observe that the curves of ^e(/i) overlap each other al- 
most completely in the Kondo regime (A ^ Tk) except 
for cases close to the mixed- valence regime {\ed\ ^ F). 
The deviation from the universal behavior in the mixed- 
valence regime (— Cd ^ F) is not surprising because of 
strong charge fluctuations in the regime. This is also in- 
dicated the phase diagram in Fig. 01 Close to the mixed- 
valence regime (— e^/F < 1), Ac/Tk becomes larger. 

More interestingly, the universal curve in Fig. ^ can fit 
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to a simple effective non-interacting model. On a non- 
interacting resonance level coupled to two superconduct- 
ing reservoirs, the spectral weights of the Andreev states 
are given by 



where 



(1 



(10) 



V{z) = z (a/i - z'^ + r) - y^e^(l-z2)-Hr2 cos2(0/2). 

(11) 

At the resonance (e^ — 0) in the limit F 3> A, the ex- 
pression is reduced to 



A 



e(h) 
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(12) 



Kondo correlated state behaves like a Fermi liquid. Nat- 
urally, if the reservoirs are normal metal, the Kondo res- 
onance can be regarded in effect as a non-interacting res- 
onance level at the Fermi energy Ep. In other words, 
many physical properties described pretty well by the ef- 
fective impurity Green's function 
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Tk/T 
z + iTk 



(13) 



with Tk playing the role of the level broadening. In 
the previous workQ with particle-hole symmetry, it was 
demonstrated that this may also work for superconduct- 
ing reservoirs. Indeed, the spectral weights in Fig. 0] fits 
very well to 



rp2 



(14) 



to be compared with Eq. H12|) . 

True bound state Finally, we address whether the 
subgap Andreev state is a true bound state. Clerk et al. 
[ssj and Sellier et al. found finite broadening of the 
Andreev states. This may come from the finite temper- 
ature effects. The NCA cannot goes down to tempera- 
tures much lower than the Kondo temperature, and they 
worked at rather high temperatures fl^ . The s pec trum 
from the NRG calculation is inherently discrete |3J|, and 
it is not easy to make a definite conclusion. However, 
as shown in Fig. |S1 the subgap states are well separated 
from the continuum parts up to temperatures as high as 
the energy of the subgap states. At temperatures higher 
than the energy of the Andreev states, it is accompa- 
nied by other small spikes. It suggests that the subgap 
Andreev states are true bound states. 

Conclusion We have studied the Kondo quantum dot 
coupled to two superconducting leads and investigated 
the subgap Andreev states using the NRG method. Con- 
trary to the recent NCA results[3, we observe An- 
dreev states both below and above the Fermi level. 
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FIG. 5; (color on-line) Raw data of the spectral weights of 
the discrete energy levels from the NRG calculation (a) at 
T = O.IA and (b) at T = 0.5. 
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